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r^wop»L^w.  ABSTRACT 

The  reversion-of-series  method  is  extended  to  the 
s  -  domain  by  utilizing  non-linear  Laplace  transforms. 
The  reversion  of  series  in  the  s  -  domain  is  applied  to  a 
non-linear  differential  equation  and  approximate  solutions 
are  obtained.   The  approximate  solution  is  modified  for 
the  case  where  the  steady  state  is  a  constant  value  by 
calculating  the  exact  steady-state  value  and  applying  it 
to  the  reversion  approximation.   The  non-linear  differen- 
tial equation  considered  is  Duffing's  equation  with  a 
damping  term  and  sinusoidal  and  constant  forcing  functions. 
The  theoretical  solutions  are  compared  to  machine  solutions 
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I.   INTRODUCTION 

Utilization  of  the  Laplace  transformation  to  obtain  an 
accurate  analytical  solution  of  a  non-linear  differential 
equation  has  been  prevented,  in  past,  by  the  inability  to 
obtain  an  appropriate  expression  for  the  Laplace  transform 
of  the  non-linear  terms.   Baycura  [Ref.  1]  developed  an 
expression  for  the  Laplace  transform  of  certain  non-linear 
terms  and  Brady  [Ref.  2]  obtained  a  formulation  for  the 
general  case.   The  non-linear  transform  is  used  in  con- 
junction with  the  reversion-of-series  method  to  obtain  an 
approximate  solution  for  a  non-linear  differential 
equation. 

The  approximate  solution  is  formulated  and  examined 
in  detail  by  comparison  to  a  machine  solution.   A  solution 
obtained  by  modifying  the  approximate  solution  is  also 
derived  and  examined  in  detail. 


II.   TRANSFORM  EXPRESSION  FOR  INTEGRAL  POWERS  OF  A  FUNCTION 


The  general  formula  for  the  Laplace  transform  of  a 
function  raised  to  an  integral  power  [Ref.  3]  is: 


^[xn(t)]  =  sn_1Xn(s) 


(1) 


where  n  is  a  positive  integer  greater  than  zero. 


III.   REVERSION  OF  SERIES 

If  a  series  is  represented  by 

2      3      4      5 
x  =  a-jY  +  a2y   +  a.^J  +   a^y   +  a^y^  +  ...  (a^O)       (2) 

the  coefficients  of  the  series 

?  3       4       5 

y  =  A,x  +  ApX   +  A^xJ  +  A^x   +  Ac.xJ+...  (3) 

are  given  by  [Ref.  4]: 

1  a2  1       2 

A1  =  —   A2  =     ^   A3  =   ^r   (2a2   -  a1a3) 

1  2        3 

Al  =   y"  (> 5&-|  a.p&^  —  a.,  a^,  —  ba^  j 

al 

A  p.  =  — q  (6a-,  apa^  +  3a-,  a.-,  +  l4a2   -  a-,  a,-  -  21a, a2  aO 
al 

Utilizing  the  relationship  expressed  in  Equation  (1)  the 
Laplace  transform  of  Equation  (2)  becomes: 

X(s)  =  a1Y(s)  +a2sY2(s)  +a3s2Y3(s)  +ai|s3Yi|(s)  +  ...      (4) 

Let: 

2  3 

b,  =  a-,    bp  =  a?s    b-  =  a~s    b^  =  a^s 

Equation  (4)  can  then  be  rewritten  as: 

X(s)  =  b1Y(s)  +  b2Y2(s)  +  b3Y3(s)  +  b^Cs)  ...  (5) 


The  Laplace  transform  of  Equation  (3)  is: 

Y(s)  =  A1X(s)  +  A2sX2(s)  +  A  s2X3(s)  +  A^s^X  (s)   ...    (6) 

Let: 

2  3 

B]  =  al    B2  =  A2S    B3  =  A3S     B4  =  A4S 

Equation  (6)  can  then  be  rewritten  as: 

Y(s)  =  B  X(s)  +  B2X2(s)  +  B  X3(s)  +  B^Cs)  ...  (7) 

The  coefficients  of  Equation  (7)  can  be  expressed  in  terms 
of  the  coefficients  of  Equation  (5): 

Bl  =  BT   B2  "  "  73  B3  =  ^5    (2b22  "  blV 

1         1 

B4   =   -^j   (5b1b2b3   -   b12b1|    -    5b23) 
bl 

Bc   =  -^77   (6b    2b0b,,    +    3b    2b02   +    l4b0      -   b    3bc    -   21b,b02b0) 

Dl 


Thus,  Y(s)  can  be  expressed  as  a  function  of  X(s)  by  the 
reversion  of  Equation  (5)  into  Equation  (7). 


IV.   GENERAL  SOLUTION  OF  A  NON-LINEAR  DIFFERENTIAL  EQUATION 

The  technique  illustrated  in  the  previous  section  will 
now  be  used  to  solve  Duffing 's  equation  with  a  damping 
term  and  initial  conditions.   The  equation  is  expressed  as: 

Hr +  att  +  woy  +  hy3  =  x(t)  (8) 

dt 

2 
where  a,  cu0  and  h  are  constants  and  x(t)  is  a  forcing  func- 
tion.  The  Laplace  transform  of  Equation  (8)  is: 

(s2  +  as  +  o)q)Y(s)  +  hs2Y3(s)  =X(s)+sy(0)  +  |yr|   +ay(0)  (9) 

Let: 

2  2 

b,  =  s   +  as  +  oon 

b3  =  hs 


(s)  =  X(s)  +  y(0)(s+a)  +  & 


dt'0 

Equation  (9)  can  then  be  written  as: 

b1Y(s)  +  b^Y3(s)  =  Z(s)  (10) 

Consider  Equation  (10)  as  series  with  b   equal  to  zero  for 
all  n  except  n  equal  to  one  and  n  equal  to  three. 
By  reversion  of  series: 

vf«1  Z(s)  hs2Z3(s)  .     3h2s1|Z5(s)  ,,,* 

7~2l      x    2^         7~2~      j.    2^4        7~2T      7    2s7  K       J 

(s    +as  +  oon)         (s    +as  +  cO  (s    +  as  +  oon) 


The  denominator  of  each  term  in  Equation  (11)  can  be 
factored.   Thus: 


2         2 
s   +  as  +  con 


s  +  ~-  +  

2        2 


a 
s  +  p- 


(12) 


Let 


n    a  ,    /  2  i.  2n 
A  =  ^  +  2"  (a  -^Wq) 


b  =  !-^(^20)2 


(13) 
(14) 


Equation  (11)  can  then  be  rewritten  as: 


Y(s)  = 


hs  Z3(s) 

U IT 

(s+A)(s+B)    (s+A)  (s+B) 


Z(s) 


3h2s4Z5(s) 
(s+A)7(s+B)7 


(15) 


Equation  (15)  represents  the  general  solution  in  the 

s  -  domain.   Inspection  of  Equation  (15)  reveals  that  the 

form  of  the  time  solution  depends  on  the  nature  of  the 

2       2 
forcing  function  and  the  relationship  between  a   and  4oon . 

Further,  if  the  forcing  function  is  Laplace  transformable 

and  if  h  is  equal  to  zero,  the  solution  is  exact. 
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V.   APPLICATION  OF  THE  GENERAL  SOLUTION 

The  general  solution  just  derived  will  now  be  applied 
for  zero  initial  conditions  with  constant  and  sinusoidal 
forcing  functions.   The  time  solutions  are  given  in  terms 
of  numbered  coefficients.   The  formula  for  each  coefficient 
is  listed  in  the  Appendix. 

A.   CONSTANT  FORCING  FUNCTION 

The  inverse  Laplace  transform  of  Equation  (15)  yields 
three  different  solutions  depending  on  the  nature  of  the 
radical  in  Equations  (13)  and  (14). 

1 .   Radical  Equal  to  Zero 

For  the  radical  to  be  equal  to  zero,  then: 

2    h    2  (16) 

and 


A  =  B  =  |  (17) 


The  general  solution  in  the  s  -  domain  is  then 


■5         2  5 
vy«o  _     c        he        3h  c  /,px 

Y(s  )  =  j  -  : 7T  +  : yh  +    •  •  •  kIo) 

s(s+A)     s(s+A)°    s(s+A) 


where  c  is  the  constant  forcing  function. 

The  inverse  Laplace  transform  of  the  first  three  terms  of 

Equation  (18)  yields: 
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y(t)  =  K]_  +  (K2+K3t)e  At  +  K^  +  (K  +Kgt+K?t2  +  .  .  .  +  Kl8t13)e  At 

(19) 

where  the  K.,  i  =  1,  2,  33  •••)  are  defined  in  the  Appendix. 

An  examination  of  Equation  (19)  indicates  that  the  initial 
conditions  are  satisfied.   The  effect  of  the  various  param- 
eters on  the  solution  can  be  seen  by  examining  the  steady- 
state  solution  which  can  be  expressed  as: 


y(t)  =  4 


to 


0  L 


1  - 


he' 


C20)3 


+  3 


he' 


C§>3 


■12 


he' 


(.§>3 


(20) 


It  should  be  noted  that  only  the  first  three  terms  of 
Equation  (20)  are  utilized  in  Equation  (19).   Since  Equation 
(20)  is  an  alternating  infinite  series,  a  condition  for 
convergence  [Ref.  5]  is  that  the  magnitude  of  each  successive 
term  must  decrease.   This  implies  that  the  applicability 

of  the  reversion  method  is  limited  in  the  steady  state  to 

2  2  3 

the  situation  where  he   is  small  compared  to  (w_)  . 

2 
Using  the  arbitrary  values  c  =  1.0,  a   =    4.0,  oon  =  4.0 

and  h  =  1.0,  the  solution  represented  by  Equation  (19) 
is  compared  in  Table  I  to  a  machine  solution  which  was 
obtained  using  Runge-Kutta  Adams  -  Moulton  with  error  check. 
As  the  results  indicate,  a  relatively  accurate  approximate 
solution  was  obtained.   In  Table  II  the  solutions  are  com- 
pared for  the  coefficient  of  the  non-linear  term  equal  to 
ten.   In  this  case  the  solution  was  a  good  approximation, 
but  not  as  accurate  as  in  the  previous  example.   Examination 
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TABLE    I 


a=4.0,    0Jq  =  4.0,    h=1.0,    c-1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0661 

0.0661 

0.0661 

0.0000 

1.0 

0.1485 

0.1485 

0.1484 

0.0001 

1.5 

0.2002 

0.2002 

0.1997 

0.0005 

2.0 

0.2271 

0.2269 

0.2259 

0.0010 

2.5 

0.2399 

0.2394 

0.2379 

0.0015 

3.0 

0.2457 

0.2447 

0.2430 

0.0017 

3.5 

0.2482 

0.2466 

0.2451 

0.0015 

4.0 

0.2492 

0.2471 

0.2458 

0.0013 

4.5 

0.2497 

0.2471 

0.2461 

0.0010 

5.0 

0.2499 

0.2467 

0.2462 

0.0005 

5;5 

0.2499 

0.2466 

0.2463 

0.0003 

6.0 

0.2500 

0.2465 

0.2463 

0.0002 

6.5 

0.2500 

0.2464 

0.2463 

0.0001 

7.0 

0.2500 

0.2463 

0.2463 

0.0000 

7.5 

0.2500 

0.2463 

0.2463 

0.0000 
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TABLE    II 


a=4.0,    co^=4.0,    h=10.0,    c=1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0661 

0.0661 

0.0660 

0.0001 

1.0 

0.1485 

0.1485 

0.1477 

0.0008 

1.5 

0.2002 

0.1997 

0.1956 

0.0031 

2.0 

0.2271 

0.2251 

0.2160 

0.0091 

2.5 

0.2399 

0.2347 

0.2221 

0.0126 

3.0 

0.2457 

0.2357 

0.2232 

0.0125 

3.5 

0.2482 

0.2327 

0.2229 

0.0098 

4.0 

0.2492 

0.2285 

0.2227 

0.0058 

4.5 

0.2497 

0.2246 

0.2225 

0.0021 

5.0 

0.2499 

0.2219 

0.2225 

-0.0006 

5.5 

0.2499 

0.2205 

0.2225 

-0.0020 

6.0 

0.2500 

0.2202 

0.2225 

-0.0023 

6.5 

0.2500 

0.2206 

0.2225 

-0.0019 

7.0 

0.2500 

0.2220 

0.2225 

-0.0005 

7.5 

0.2500 

0.2233 

0.2225 

0.0008 

8.0 

0.2500 

0.2246 

0.2225 

0.0021 

9.0 

0.2500 

0.2268 

0.2225 

0.0043 

15.0 

0.2500 

0.2292 

0.2225 

0.0067 
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of  Equation  (20)  reveals  that  for  h  equal  to  ten,  successive 
terms  decrease  in  magnitude,  but  not  as  rapidly  as  in  the 
case  h  was  equal  to  one.   If  four  terms  of  Equation  (20) 
had  been  utilized  in  the  solution,  a  steady-state  value  of 
0.2178  would  have  been  obtained  resulting  in  a  better 
approximation.   Table  III  represents  a  comparison  of  the 
solutions  for  the  forcing  function  equal  to  ten.   In  this 
case  the  solution  limits  were  exceeded  and  a  large  positive 
number  was  obtained  for  the  steady  state.   If  an  additional 
term  of  Equation  (20)  had  been  included  in  the  solution,  a 
larger  negative  value  would  have  been  obtained  for  the 
steady  state.   Good  results  were  realized  however,  for  the 
initial  portion  of  the  transient  solution,  indicating  the 
applicability  of  the  reversion  method  in  this  region.   The 
linear  solution  represented  by  the  first  two  terms  of 
Equation  (19)  was  also  included  in  the  tables  to  illustrate 
a  less  accurate  but  simplified  approximate  solution. 

2 .   Radical  Greater  than  Zero 

2       2 
The  relationship  between  a   and  4oon  is  expressed  as: 

a2>4WQ  .  (21) 

Equation  (15)  is  then  expressed  as: 

Y(S)  = s jjsi! — +  — 3h!ci —  (22) 

s(s+A)(s+B)    s(s+A)  (s+B)H   s ( s+A) ' ( s+B) ' 

The  inverse  Laplace  transform  of  the  first  two  terms  of 
Equation  (22)  yields: 
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TABLE    III 


a=4.0,    o3q=4.0,    h=1.0,    c  =  10.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.2 

0.1539 

0.1539 

0.1539 

0.0000 

0.4 

0.4780 

0.4780 

0.4770 

0.0010 

0.6 

0.8434 

0.8433 

0.8361 

0.0072 

0.8 

1.1877 

1.1867 

1.1628 

0.0239 

1.0 

1.4850 

1.4807 

1.4085 

0.0722 

1.5 

2.0020 

1.9560 

1.6365 

0.3195 

2.0 

2.2710 

2.0727 

1.5880 

0.4847 

2.5 

2.3990 

1.8907 

1.5006 

0.3901 

3.0 

2.4570 

1.5230 

1.5515 

-0.0285 

3.5 

2.4818 

1.1488 

1.5567 

-0.4079 

4.0 

2.4925 

0.9814 

1.5575 

-0.5761 

4.5 

2.4969 

1.2081 

1.5569 

-0.3488 

5.0 

2.4988 

1.9344 

1.5567 

0.3777 

6.0 

2.4998 

4.7740 

1.5568 

3.2172 

7.0 

2.5000 

8.5235 

1.5568 

6.9967 

8.0 

2.5000 

11.9169 

1.5568 

10.3601 

9.0 

2.5000 

14.3046 

1.5568 

12.7478 

10.0 

2.5000 

15.6965 

1.5568 

14.1397 

20.0 

2.5000 

16.9043 

1.5568 

15.3475 
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yCt)  =  K2Q  +K21e  At  +  K22e~Bt  +  K?3  +  (K^  +  YL      t   +  K^t2 +K2?t3)e  At 
+  (K2g  +  K29t  +K3Qt2  +K31t3)e"Bt  .  (23) 


The  inverse  Laplace  transform  of  additional  terms  of 
Equation  (22)  results  in  a  steady-state  solution  given  by 
Equation  (20).   The  region  in  which  the  solution  can  be 
applied  is  then  the  same  as  discussed  in  the  previous 
section.   Using  the  values  c  =  1.0,  a   =  6.0,  con  =  4.0  and 
h  =  1.0,  the  solution  represented  by  Equation  (23)  and  the 
solution  obtained  considering  the  first  three  terms  of 
Equation  (23)  are  compared  to  a  machine  solution  in  Table 
IV.   The  small  errors  listed  indicate  that  a  highly  accurate 
approximation  was  obtained.   Table  V  illustrates  the  effect 
of  increasing  the  coefficient  of  the  non-linear  term  to  ten. 
A  less  accurate  approximation  was  obtained.   The  effect  of 
increasing  the  magnitude  of  the  forcing  function  was  not 
illustrated;  however  for  the  parameter  values  used  in 
Table  IV  if  c  had  been  increased  to  ten,  a  steady-state 
value  of  -1.4063  would  have  been  obtained  indicating  that 
the  solution  limits  would  have  been  exceeded.   The  linear 
solution  is  seen  to  be  a  less  accurate  approximation. 

3 .   Radical  Imaginary 

2  2 

The    relationship   between    a      and    4ton    is: 

4oo2    >    a2    .  (24) 
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TABLE    IV 


a=6.0,    0Jq  =  4.0,    h=1.0,    c=1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0533 

0.0533 

0.0533 

0.0000 

1.0 

0.1139 

0.1139 

0.1138 

0.0001 

1.5 

0.1570 

0.1569 

0.1568 

0.0001 

2.0 

0.1865 

0.1864 

0.1860 

0.0004 

2.5 

0.2066 

0.2064 

0.2058 

0.0006 

3.0 

0.2204 

0.2200 

0.2191 

0.0009 

3.5 

0.2298 

0.2292 

0.2281 

0.0011 

4.0 

0.2362 

0.2353 

0.2341 

0.0012 

5.0 

0.2436 

0.2420 

0.2408 

0.0012 

6.0 

0.2470 

0.2448 

0.2438 

0.0010 

7.0 

0.2486 

0.2459 

0.2452 

0.0007 

8.0 

0.2494 

0.2462 

0.2459 

0.0003 

9.0 

0.2497 

0.2463 

0.2460 

0.0003 

10.0 

0.2499 

0.2463 

0.2462 

0.0001 

11.0 

0.2499 

0.2462 

0.2462 

0.0000 

12.0 

0.2500 

0.2462 

0.2463 

-0.0001 

13.0 

0.2500 

0.2461 

0.2463 

-0.0002 
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TABLE   V 


a=6.0,    0^=4.0,    h=10.0,    c=1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0533 

0.0533 

0.0533 

0.0000 

1.0 

0.1139 

0.1139 

0.1136 

0.0003 

1.5 

0.1570 

0.1568 

0.1553 

0.0015 

2.0 

0.1865 

0.1857 

0.1821 

0.0036 

2.5 

0.2066 

0.2047 

0.1986 

0.0061 

3.0 

0.2204 

0.2165 

0.2086 

0.0079 

3.5 

0.2298 

0.2234 

0.2145 

0.0089 

4.0 

0.2362 

0.2269 

0.2179 

0.0090 

4.5 

0.2406 

0.2281 

0.2196 

0.0085 

5.0 

0.2436 

0.2278 

0.2210 

0.0068 

5.5 

0.2465 

0.2266 

0.2216 

0.0050 

6.0 

0.2470 

0.2250 

0.2220 

0.0030 

6.5 

0.2480 

0.2231 

0.2222 

0.0009 

7.0 

0.2486 

0.2213 

0.2223 

-0.0010 

8.0 

0.2494 

0.2181 

0.2224 

-0.0033 

9.0 

0.2497 

0.2157 

0.2225 

-0.0068 

10.0 

0.2499 

0.2139 

0.2225 

-0.0086 

12.0 

0.2500 

0.2120 

0.2225 

-0.0109 

15.0 

0.2500 

0.2110 

0.2225 

-0.0115 
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Let 


\   (a2-^2)  2  =  jD 


Then: 


A  =  |  +  jD  (25) 


and 


B  =  |  -  JD  .  (26) 

Equation  (15)  is  then  expressed  as: 

3  oi2>3 


he  3h  c 

s(s+A)(s+B)    s(s+AT(s+Br    s  (  s+A)  7(s+B)  7 


Y(s)  =  * ^ jj-  +  ^V *  +  ... 


(27) 


The  inverse  Laplace  transform  of  the  first  two  terms  of 
Equation  (27)  yields: 

y(t)  =K40+K4le-At  +K42e-Bt  +  K43+  (K^  +  K^t  +  K^t^K^t3)  e~At 

+  (K4g  +K49t  +K5Qt2  +K51t3)e~Bt.  (28) 

Since  A  and  B  are  complex,  the  coefficients  of  the  time- 
dependent  portion  of  Equation  (28)  are  complex.   Equation 
(28)  can  then  be  rewritten  as: 

y(t)  =  K4q  +  (M1+jM2)e"At  +  (I^-jM^e-^  +  K^ 

?  3        -  At" 

+  [(M3+jM4)   + (M5+jM6)t  +  (M7+jMg)t"  +  (Mg+jM10)tJ]e 

+[(M3-jM4)   +(M5-jM6)t  +  (M?-jM8)t2   +(Mg-jM10)t3]e"Bt. 

(29) 
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Combining  terms  Equation  (29)  becomes: 

-— t  -— t 

y(t)  =  K^q  +  2M1e  2  cos(Dt)  +  2M2e  2  sin(Dt) 

+    (2M0+2Mc:t  +  2Mr7t2  +  2Mnt3)e    2    cos(Dt) 

--t 

+  (2Mi4  +  2M6t  +  2Mgt2  +  2M10t3)e  2  sin(Dt)  .  (30) 


In  this  case  as  in  the  previous  two  cases  considered,  the 
steady-state  solution  is  given  by  Equation  (20).   In  Table 
VI,  the  reversion  solution  is  compared  to  a  machine  solution 
for  a  =  1.0,  a)2  =  4.0,  h  =  1.0  and  c  =  1.0.   In  Table  VII, 
the  solutions  are  compared  for  the  non-linear  coefficient 
equal  to  ten.   A  comparison  of  the  results  listed  in  Tables 

VI  and  VII  with  the  results  obtained  in  the  previous  cases 

2 
for  the  same  values  of  con,  h  and  c,  indicates  that  the 

reversion  method  yields  a  less  accurate  transient  solution 

when  the  damping  coefficient  is  small.   In  this  case  the 

linear  solution  is  seen  to  be  a  better  approximation  over  a 

portion  of  the  transient  solution. 

B.   SINUSOIDAL  FORCING  FUNCTION 

Substitution  of  a  constant-amplitude  sinusoidal  forcing 
function,  c  sin  (cot),  into  the  first  two  terms  of  Equation 
(15)  and  utilization  of  the  trigonometric  identity  [Ref.  61 

■DO  ■» 

sin  (cot)  =  -jsin(cot)  -  £-sin(  3^t )  (31) 

yields : 
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TABLE    VI 


a=1.0,    0)^  =  4.0,    h=1.0,    c  =  1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0982 

0.0982 

0.0982 

0.0000 

1.0 

0.2677 

0.2677 

0.2672 

0.0005 

1.5 

0.3576 

0.3574 

0.3534 

0.0040 

2.0 

0.3343 

0.3330 

0.3235 

0.0095 

2.5 

0.2591 

0.2550 

0.2463 

0.0087 

3.0 

0.2069 

0.1985 

0.1994 

-0.0009 

3.5 

0.2064 

0.1951 

0.2060 

-0.0109 

4.0 

0.2377 

0.2279 

0.2402 

-0.0123 

4.5 

0.2655 

0.2620 

0.2655 

-0.0035 

5.0 

0.2712 

0.2743 

0.2663 

0.0080 

5.5 

0.2593 

0.2642 

0.2516 

0.0126 

6.0 

0.2453 

0.2453 

0.2387 

0.0066 

6.5 

0.2403 

0.2324 

0.2367 

-0.0043 

7.0 

0.2442 

0.2320 

0.2426 

-0.0106 

8.0 

0.2542 

0.2510 

0.2508 

0.0002 

9.0 

0.2503 

0.2523 

0.2454 

0.0069 

10.0 

0.2483 

0.2409 

0.2450 

-0.0041 

50.0 

0.2500 

0.2461 

0.2463 

-0.0002 
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TABLE   VII 


a=1.0,    to^=4.0,    h=10.0,    c-1,0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.0000 

0.5 

0.0982 

0.0982 

0.0982 

0.0000 

1.0 

0.2677 

0.2675 

0.2626 

0.0049 

1.5 

0.3576 

0.3554 

0.3188 

0.0366 

2.0 

0.3343 

0.3211 

0.2507 

0.0704 

2.5 

0.2591 

0.2173 

0.1799 

0.0374 

3.0 

0.2069 

0.1230 

0.1761 

-0.0531 

3.5 

0.2064 

0.0934 

0.2164 

-0.1230 

4.0 

0.2377 

0.1400 

0.2463 

-0.1063 

4.5 

0.2655 

0.2298 

0.2410 

-0.0112 

5.0 

0.2712 

0.3026 

0.2195 

0.0831 

5.5 

0.2593 

0.3085 

0.2093 

0.0992 

6.0 

0.2453 

0.2451 

0.2156 

0.0295 

6.5 

0.2403 

0.1619 

0.2260 

-0.0641 

7.0 

0.2442 

0.1222 

0.2289 

-0.1067 

8.0 

0.2542 

0.2225 

0.2198 

0.0027 

9.0 

0.2503 

0.2703 

0.2222 

0.0481 

10.0 

0.2483 

0.1745 

0.2234 

-0.0489 

50.0 

0.2500 

0.2110 

0.2225 

-0.0115 
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3   3 

„(  ■y     _    uc _  4 , 

(s2+to2)(s+A)(s+B)    (s2+w2)(s+A)4(s+B)4 


3   3 

fhc  a) 


(s2  +  9oo2)(s+A)4(s+B)i| 


(32) 


As  in  the  case  of  the  constant  forcing  function  three  solu- 
tions are  considered  depending  on  the  nature  of  the  radical 
in  Equations  (13)  and  (14). 

1 .   Radical  Greater  than  Zero 

The  inverse  Laplace  transform  of  Equation  (32)  is: 

y(t)  =  M11cos(o)t)  +  M12sin(u)t)  +  N1e~At  +  N2e"Bt+  M-^cosCwt)  + 

Ml4sin(o3t)  +  (N_  +  N^t  +  N5t2  +  N6t3)e~At  + 

(N„  +  Ngt  +  Ngt2  +  N10t3)e"Bt  +  M   cos(3wt)  +  Ml6sin(3wt)  + 

(N1X  +  N12t  +  N13t2  +Nx^t3)e"At+  (N15+N^t  +  N17t2  +  Nl8t3)e"Bt. 

(33) 
The  solutions  represented  by  Equation  (33)  and  the  first 

four  terms  of  Equation  (33)  were  tested  using  a  =  8.0, 

2 
(j0n=9.0,  c  =  1.0,  h=20.0  and  to=1.0.   The  solutions  are  compared 

to  a  machine  solution  in  Table  VIII.  In  Table  IX  the  solu- 
tions are  compared  for  the  non-linear  coefficient  increased 
to  fifty.   Table  X  is  a  comparison  of  the  solutons  for  the 
frequency  of  the  forcing  function  increased  to  five  and  the 
non-linear  coefficient  equal  to  fifty.   Inspection  of  the 
tables  indicates  that  as  the  non-linear  coefficient  is 


24 


TABLE   VIII 
ot=8.0,    a)^=9.0}    h=20.0,    c=1.0,   oo  =  1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.00000 

0.00000 

0.00000 

0.00000 

0.5 

0.00885 

0.00885 

0.00885 

0.00000 

1.0 

0.03603 

0.03603 

0.03603 

0.00000 

1.5 

0.06667 

0.06666 

0.06661 

0.00005 

2.0 

0.08728 

0.08726 

0.08704 

0.00022 

2.5 

0.08973 

O.O8963 

0.08923 

0.00040 

3.0 

0.07184 

0.07160 

0.07128 

0.00032 

3.5 

0.03719 

0.03679 

0.03674 

0.00005 

4.0 

-0.00615 

-0.00667 

-0.00644 

0.00023 

4.5 

-0.04777 

-0.04833 

-0.04789 

0.00044 

5.0 

-0.07759 

-0.07810 

-0.07753 

0.00057 

5.5 

-0.08835 

-0.08873 

-0.08803 

0.00070 

6.0 

-0.07745 

-0.07762 

-0.07695 

0.00067 

6.5 

-0.04758 

-0.04750 

-0.04714 

0.00036 

7.0 

-0.00605 

-0.00576 

-0.00579 

-0.00003 

23.0 

-0.01959 

-0.01998 

-0.01991 

0.00007 

24.0 

-0.08311 

-0.08352 

-0.08298 

0.00054 

25.0 

-0.07022 

-0.07029 

-0.06971 

0.00058 

26.0 

0.00723 

0.00758 

0.00745 

0.00013 

27.0 

0.07803 

0.07847 

0.07798 

0.00049 

28.0 

0.07709 

0.07722 

0.07658 

0.00064 

29.0 

0.00528 

0.00497 

0.00502 

-0.00005 
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TABLE    IX 


a=8.0,    a)^=9.0,    h=50.0,    c  =  1.0,w=1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.00000 

0.00000 

0.00000 

0.00000 

0.5 

0.00885 

0.00885 

0.00885 

0.00000 

1.0 

0.03603 

0.03603 

0.03601 

0.00002 

1.5 

0.06667 

0.06665 

0.06637 

0.00028 

2.0 

0.08728 

0.08715 

0.08606 

0.00109 

2.5 

0.08973 

0.08922 

0.08732 

0.00190 

3.0 

0.07184 

0.07064 

0.06905 

0.00059 

3.5 

0.03719 

0.03519 

0.03510 

0.00009 

4.0 

-0.00615 

-0.00874 

-0.00731 

0.00143 

4.5 

-0.04777 

-0.05056 

-0.04831 

0.00225 

5.0 

-0.07759 

-0.08015 

-0.07730 

0.00085 

5.5 

-0.08835 

-0.09025 

-0.08677 

0.00348 

6.0 

-0.07745 

-0.07829 

-0.07508 

0.00321 

6.5 

-0.04758 

-0.04717 

-0.04544 

0.00173 

7.0 

-0.00605 

-0.00458 

-0.00479 

-0.00021 

22,8 

-0.00207 

-0.00375 

-0.00326 

0.00049 

24.0 

-0.08311 

-0.08517 

-0.08243 

0.00274 

25.0 

-0.07022 

-0.07054 

-0.06780 

0.00274 

26.0 

0.00723 

0.00899 

0.00831 

0.00068 

27.0 

0.07803 

0.08020 

0.07769 

0.00251 

28.0 

0.07709 

0.07774 

0.07465 

0.00309 

29.0 

0.00528 

0.00373 

0.00402 

-0.00029 
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TABLE    X 


a=8.0,    w^=9.0,    h=50.0,    c=1.0,    00=5.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.00000 

0.00000 

0.00000 

0.00000 

0.2 

0.0043^ 

.   0.00434 

0.00434 

0.00000 

0.4 

0.02066 

0.02066 

0.02066 

0.00000 

0.6 

0.03549 

0.03549 

0.03549 

0.00000 

0.8 

0.03246 

0.03246 

0.03242 

0.00004 

1.0 

0.01123 

0.01121 

0.01115 

0.00006 

1.2 

-0.01136 

-0.01140 

-0.01144 

-0.00004 

1.4 

-0.01662 

-0.01672 

-0.01670 

-0.00002 

1.6 

-0.00136 

-0.00151 

-0.00141 

0.00010 

1.8 

0.01916 

0.01895 

0.01913 

-0.00018 

2.0 

0.02512 

0.02487 

0.02511 

-0.00024 

2.2 

0.01032 

0.01005 

0.01029 

-0.00024 

2.4 

-0.01220 

-0.01248 

-0.01224 

0.00024 

2.6 

-0.02214 

-0.02243 

-0.02218 

0.00025 

2.8 

-0.01069 

-0.01100 

-0.01072 

0.00028 

3.0 

0.01137 

0.01106 

0.01137 

-0.00031 

6.0 

0.00520 

0.00515 

0.00519 

-0.00004 

6.4 

-0.02273 

-0.02276 

-0.02275 

0.00001 

6.8 

0.01373 

0.01370 

0.01374 

0.00004 

7.2 

0.01131 

0.01129 

0.01131 

0.00003 

7.6 

-0.02314 

-0.02314 

-0.02316 

0.00002 
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increased,  the  accuracy  of  the  solution  decreases  and  as 
the  frequency  of  the  forcing  function  is  increased,  a 
more  accurate  approximation  is  obtained.   The  effect  of  the 
magnitude  of  the  amplitude  of  the  forcing  function  was  not 
illustrated;  however  from  an  examination  of  the  coefficients 
of  Equation  (33)  it  can  be  seen  that  as  c  is  increased  it 
would  have  a  more  detrimental  effect  on  the  accuracy  of 
the  solution  than  corresponding  increases  of  the  non-linear 
coefficient.   An  examination  of  the  solution  for  the  non- 
linear coefficient  equal  to  zero,  represented  by  the  first 
four  terms  of  Equation  (33)>  reveals  that  it  is  overall  as 
accurate  a  solution  as  the  solution  utilizing  all  the  terms 
of  Equation  (33) • 

2 .   Radical  Equal  to  Zero 
Equation  (32)  becomes: 


11  11 

tfhc    oo  -jThc    0) 

y  (     S      ,    coc 4 4 (  r>u\ 

^s '              o      o                o             o      o                S"  oo                p    "    ^  -3    ' 

(s^+oj    Ks+AT         (s^+u)    )(s+A)°  (s^  +  9oj    )(s  +  A)b 


The    inverse   Laplace    transform   of   Equation    (34)    yields: 
y(t)    =   M17cos(a)t)  +Ml8sin(ajt)   +  (N19+  N20t)e"At+   M19cos(ajt)    + 
M20sinCwt)  +  CN21+N22t+N2    t2+N24t3+N25t^+N26t5+N2    t6+N28t7)e~At+ 
M21cos(3ut)+M22sin(3wt)  +  (N      +N      t+N      t2+N      t3+N„t      + 
NOJlt5  +  Nocrt6  +N0.t7)e"At    .  (35) 


34  35  3d 
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Table  XI  is  a  comparison  of  the  solutions  for  the  parameter 
values  a=6.0,  con=9.0,  h=20.03  c=1.0  and  oj=1.0.   A  comparison 
of  errors  listed  in  Table  XI  with  the  errors  listed  in 
Table  VIII  indicates  that  a  decrease  in  the  magnitude  of 
the  damping  coefficient  results  in  a  less  accurate  approxi- 
mate solution.   The  linear  solution  represented  by  the 
first  three  terms  of  Equation  (35)  is  again  seen  to  be 
overall  as  accurate  an  approximation  when  compared  to  the 
solution  represented  by  Equation  (35). 
3 .   Radical  Imaginary 

Utilizing  the  relationships  expressed  in  Equations 
(25)  and  (26)  and  combining  terms  with  complex  conjugate 
coefficients j  the  inverse  Laplace  transform  of  Equation  (32) 

becomes : 

a,  au. 

y(t)  =  M2-cos(u>t)  +M  i|sin((ut)  +N^e    cos  (Dt  )+N^ge   sin(Dt)  + 

-St 

M2C-cos(tot)+M2gSin(ajt  )  +  (N^q+Ni40t+NI|1t2+Ni|2t3)e    cos(Dt)  + 

-^t 
(Nij^+Ni4ijt+Ni|C-t2+N1|6t3)e  2  sin(Dt)+M27cos(3wt)  + 

-ft 

M2gsin(3wt  )  +  (Ni|7+Ni|gt+N1|9t2+Nt-0t3)e  2  cos(Dt)  + 

-^t 
(N  +N  t+N  t2+N5i|t3)e  2  sin(Dt)  .  (36) 

Using  the  values  a=1.0,  u)?  =  9.0,  h=20.0,  c-1.0  and  uj=1.0, 
the  solutions  are  compared  in  Table  XII.   A  comparison  of 
Table  XII  with  Tables  VIII  and  XI  further  illustrates  the 
critical  effect  of  the  damping  coefficient  on  the  accuracy 
of  the  solution.  In  the  three  tables  the  parameter  values  are 
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TABLE    XI 


a=6.0,    a)^  =  9.0,    h=20.0,    c=1.0,    w=1.0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.00000 

0.00000 

0.00000 

0.00000 

0.5 

0.01024 

0.01024 

0.01025 

0.00001 

1.0 

0.04287 

0.04287 

0.04285 

0.00002 

1.5 

0.07789 

0.07788 

0.07766 

0.00022 

2.0 

0.09836 

0.09825 

0.09746 

0.00079 

2.5 

0.09612 

0.09572 

0.09449 

0.00123 

3.0 

0.07073 

0.06984 

0.06912 

0.00072 

3.5 

0.02814 

0.02675 

0.02714 

-0.00039 

4.0 

-0.02132 

-0.02290 

-0.02171 

0.00119 

4.5 

-0.06555 

-0.06695 

-0.06560 

0.00135 

5.0 

-0.09373 

-0.09470 

-0.09320 

0.00140 

5.5 

-0.09896 

-0.09932 

-0.09759 

0.00173 

6.0 

-0.07996 

-0.07958 

-0.07827 

0.00131 

6.5 

-0.04139 

-0.04030 

-0.04020 

0.00010 

7.0 

0.00732 

O.OO88O 

0.00783 

0.00097 

21.0 

0.09980 

0.10042 

0.09875 

0.00167 

22.0 

0.05929 

0.05847 

0.05782 

0.00065 

23.0 

-0.03573 

-0.03723 

-0.03600 

0.00123 

24.0 

-0.09790 

-0.09869 

-0.09711 

0.00158 

25.0 

-0.07006 

-0.06945 

-0.06846 

0.00099 

26.0 

0.02219 

0.02370 

0.02256 

0.00114 

27.0 

0.09404 

0.09498 

0.09349 

0.00149 
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TABLE    XII 
01=1.0,    u>q  =  9.0,    h=20.0,    c=1.0,    oo=l .  0 


Time 

Linear 

Reversion 

Machine 

Error 

0.0 

0.00000 

0.00000 

0.00000 

0.00000 

0.5 

0.01635 

0.01635 

0.01635 

0.00000 

1.0 

0.08176 

0.08174 

0.08161 

0.00013 

1.5 

0.13744 

0.13703 

0.13485 

0.00218 

2.0 

0.12875 

0.12637 

0.12134 

0.00503 

2.5 

0.07792 

0.07206 

0.07300 

0.00094 

3.0 

0.02512 

0.01893 

0.02743 

-0.00850 

3.5 

-0.02494 

-0.02464 

-0.02176 

0.00288 

4.0 

-0.07801 

-0.07121 

-0.07877 

-0.00756 

4.5 

-0.11866 

-0.11433 

-O.H872 

-0.00439 

5.0 

-0,12570 

-0.12773 

-0.12008 

0.00765 

5.5 

-0.09724 

-0.09750 

-0.09118 

0.00632 

6.0 

-0.04708 

-0.04019 

-0.04717 

-0.00698 

6.5 

0.01145 

0.01727 

0.00821 

0.00906 

7.0 

0.06800 

0.06440 

0.06752 

-0.00312 

22.0 

0.01429 

0.01258 

0.01643 

-0.00385 

23.0 

-0.09595 

-0.09420 

-0.09641 

-0.00221 

24.0 

-0.11798 

-0.11444 

-0.11215 

0.00229 

25.0 

-0.03154 

-0.02945 

-0.03276 

-0.00331 

26.0 

0.08390 

0.08262 

0.08404 

-0.00142 

27.0 

0.12220 

0.11872 

0.11709 

0.00163 

28.0 

0.04815 

0.04568 

0.04806 

-0.00238 
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equal  except  for  the  magnitude  of  the  damping  coefficient. 
As   a  was  decreased  from  eight  in  Table  v"III  to  six  in 
Table  XI  and  finally  to  one  in  Table  XII',   the  magnitude 
of  the  error  increased  correspondingly.   The  linear  solu- 
tion illustrated  in  Table  XII  is  seen  to  be  an  overall 
better  approximation  than  the  reversion  solution. 
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VI.   MODIFIED  SOLUTION 

The  reversion  method  previously  examined  was  limited  in 
application  to  a  region  in  which  the  solution  did  not 
diverge.   The  reversion  method  can  be  made  to  always  con- 
verge by  modifying  the  solution  to  be  correct  in  the  steady 
state  for  any  parameter  values.   The  steady  state  for  a 
constant  forcing  function  can  be  obtained  by  solving  the 
cubic  equation: 

w^y  +  hy3  =  c  .  (37) 

Consider   a   truncated   form  of  Equation    (23): 

y(t)    =    K2Q  +K21e~At  +K22e"Bt  +K23    +  (K2l|  +  K^t  )e"At    + 

<K28  +  K29t)e-Bt    •  (38) 

Since  the  coefficients  K?~  to  K?q  have  the  common  factor 
he  ,  see  Appendix,  Equation  (38)  can  be  rewritten  as: 

3 
y(t)  =  K2Q  +  K21e"At  +  K22e~Bt  -  h°    k    +  hc3(K24  +  K2  t)e~At  + 

hc3(K28  +K29t)e"Bt  .  (39) 

Let  SS  be  equal  to  the  steady  state  obtained  from  the  solu- 
tion of  Equation  (37).   If  the  steady-state  portion  of 
Equation  (39)  is  set  equal  to  SS,  then: 

ss  =  K20  ~  TTT^  '  (40) 

C»0) 
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Equation  (40)  can  be  rewritten  as: 

he3  =  (o)q)4(K20  -  SS)  (4l) 

Substitution  of  Equation  (41)  into  Equation  (39)  yields: 
y(t)  =  SS+K21e  RZ  +   K^e  OL  +  (ug)  H  (K2Q-SS)  (K2l|  +  K^t  )e  At  + 

Pi]  '       '      -Rt" 

(o3q)4(K20  -SS)(K2g  +  K29t)e     .  (42) 

Equation  (42)  is  the  modified  solution.   Using  the  values 
a=6.0,  oo  =4.0,  h=1.0  and  c=1.0,  the  modified  solution  is 
compared  to  a  machine  solution  in  Table  XIII.   The  results 
indicate  a  relatively  accurate  approximate  solution  was 
obtained.   In  Tables  XIV  and  XV  the  solutions  are  compared 
for  the  non-linear  coefficient  respectively  equal  to  ten 
and  fifty.   The  tables  reveal  that  the  error  increased  as  h 
was  increased;  however,  the  modified  solution  is  seen  to 
be  still  usable  as  an  approximation  for  the  non-linear 
coefficient  equal  to  fifty.   In  Table  XVI  the  solutions  are 
compared  for  the  forcing  function  equal  to  ten.   The  magni- 
tude of  the  forcing  function  is  seen  to  be  a  more  critical 
factor  on  the  accuracy  of  the  approximation  than  correspond- 
ing increases  of  the  non-linear  coefficient.   A  comparison 
of  the  modified  solution  with  the  reversion  solution  for 
similar  parameter  values  indicates  that  the  modified  solution 
is  a  more  accurate  approximation  for  large  values  of  h  and  c. 
For  small  values  of  h  and  c,  the  modified  and  reversions 
solutions  are  of  the  same  relative  accuracy,  as  revealed  by 
a  comparison  of  Table  XIII  with  Table  IV. 
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TABLE    XIII 


a=6.0,    03^=4.0,    h=1.0,    c  =  1.0 


rime 

Modified 

Machine 

0.0 

0.0000 

0.0000 

0.5 

0.0532 

0.0533 

1.0 

0.1133 

0.1138 

1.5 

0.1559 

0.1568 

2.0 

0.1849 

0.1860 

2.5 

0.2046 

0.2058 

3.0 

0.2180 

0.2191 

3.5 

0.2271 

0.2281 

4.0 

0.2333 

0.2341 

4.5 

0.2374 

0.2381 

5.0 

0.2403 

0.2408 

5.5 

0.2422 

0.2426 

6.0 

0.2435 

0.2438 

6.5 

0.2444 

0.2446 

7.0 

0.2450 

0.2452 

8.0 

0.2457 

0.2459 

9.0 

0.2460 

0.2461 

10.0 

0.2461 

0.2462 

12.0 

0.2462 

0.2463 

14.0 

0.2463 

0.2463 

Error 


0.0000 
■0.0001 
■0.0005 
■0.0011 
-0.0011 
■0.0012 
-0.0011 
-0.0010 
-0.U008 
-0.0007 
-0.0005 
-0.0004 
-0.0003 
-0.0002 
-0.0002 
-0.0002 
-0.0001 
-0.0001 
-0.0001 

0.0000 
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TABLE    XIV 


a=6.0,    cOq  =  4.0,    h=10.0,    c=1.0 


Time 

Modified 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.5 

0.0522 

0.0533 

-0.0011 

1.0 

0.1098 

0.1136 

-0.0038 

1.5 

0.1493 

0.1553 

-0.0060 

2.0 

0.1751 

0.1821 

-0.0070 

2.5 

0.1919 

0.1986 

-0.0067 

3.0 

0.2028 

0.2086 

-0.0058 

3.5 

0.2099 

0.2145 

-0.0046 

4.0 

0.2145 

0.2179 

-0.0034 

4.5 

0.2174 

0.2196 

-0.0022 

5.0 

0.2193 

0.2210 

-0.0017 

5.5 

0.2205 

0.2216 

-0.0011 

6.0 

0.2212 

0.2220 

-0.0008 

6.5 

0.2217 

0.2222 

-0.0005 

7.0 

0.2220 

0.2223 

-0.0003 

7.5 

0.2222 

0.2224 

-0.0002 

8.0 

0.2223 

0.2224 

-0.0001 

9.0 

0.2224 

0.2225 

-0.0001 

10.0 

0.2225 

0.2225 

0.0000 
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TABLE   XV 


a=6.0i    co^  =  4.0,    h=50.0,    c=1.0 


Time 

Modified 

Machine 

Error 

0.0 

0.0000 

0.0000 

0.0000 

0.5 

0.0503 

0.0533 

-0.0030 

1.0 

0.1034 

0.1123 

-0.0089 

1.5 

0.1370 

0.1489 

-0.0119 

2.0 

0.1570 

0.1671 

-0.0101 

2.5 

0.1685 

0.1747 

-0.0062 

3.0 

0.1749 

0.1774 

-0.0025 

3.5 

0.1782 

0.1783 

-0.0001 

4.0 

0.1799 

0.1786 

0.0013 

4.5 

0.1805 

0.1787 

0.0018 

5.0 

0.1806 

0.1787 

0.0019 

5.5 

0.1805 

0.1787 

0.0018 

6.0 

0.1802 

0.1787 

0.0015 

8.0 

0.1793 

0.1787 

0.0006 

10.0 

0.1789 

0.1787 

0.0002 

12.0 

0.1787 

0.1787 

0.0000 
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TABLE   XVI 


a=6.0,    u)?  =  4.0,    h=1.0,    c  =  10.0 


Time 

Modified 

Machine 

0.0 

0.0000 

0.0000 

0.5 

0.4929 

0.5326 

1.0 

1.0000 

1.1077 

1.5 

1.3061 

1.4187 

2.0 

1.4749 

1.5294 

2.5 

1.5619 

1.5545 

3.0 

1.6019 

1.5574 

3.5 

1.6160 

1.5571 

4.0 

1.6166 

1.5568 

4.5 

1.6093 

1.5568 

5.0 

1.6027 

1.5568 

6.0 

1.5866 

1.5568 

7.0 

1.5746 

1.5568 

8.0 

1.5669 

1.5568 

10.0 

1.5598 

1.5568 

12.0 

1.5576 

1.5568 

14.0 

1.5570 

1.5568 

16.0 

1.5568 

1.5568 

Error 


0.0000 
•0.0397 
■0.1077 
■0.1126 

-0.0545 
0.0074 
0.0445 
0.0589 
0.0598 
0.0525 
0.0459 
0.0298 
0.0178 
0.0101 
0.0030 
0.0008 
0.0002 
0.0000 


VII.   CONCLUSIONS 

A  method  for  obtaining  an  approximate  solution  for  a 
specific  non-linear  differential  equation  has  been  derived 
utilizing  reversion-of-series  and  non-linear  transform 
techniques.   The  usefulness  of  the  solution  was  found  to 
be  dependent  upon  the  nature  of  the  forcing  function  and 
the  parameter  values  specified.   Also,  it  was  found  that 
useful  approximation  could  be  obtained  for  particular  par- 
ameter values  by  simply  considering  as  a  solution  that 
portion  of  the  reversion  approximation  which  was  independent 
of  the  non-linear  coefficient.   For  the  case  of  the  constant 
forcing  function  the  linear  approximation  was  found  to  be 
less  accurate  than  the  reversion  solution.   In  the  case 
of  the  sinusoidal  forcing  function  the  linear  approximation 
was  found  to  be  as  accurate  as  the  reversion  solution  when 
the  damping  was  significant  and  a  more  accurate  approxima- 
tion when  the  damping  coefficient  was  relatively  small. 

For  a  constant  forcing  function  it  was  found  that  the 
steady  state  could  be  specified  and  utilized  to  modify  the 
reversion  solution.   The  modified  solution  was  found  to  be 
useful  as  an  approximation  over  a  relatively  wide  range  of 
parameter  values  and  gave  results  that  were  superior  to  the 
linear  and  reversion  approximations. 

The  reversion  method  was  used  to  solve  a  particular  non- 
linear differential  equation;  however,  the  development  of  the 
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solution  indicates  the  possibility  of  applying  this  method 
to  solve  non-linear  differential  equations  similar  in  form 
to  the  equation  considered. 
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K^ 


APPENDIX 


COEFFICIENT  FORMULAS 


K„  =  -  £ 


K3  = 


A 


K^  » 


3      2  5 
he-3    3h  cD 


K5  =  -  K4 


K   = 
n 


K   = 
n 


1 


(n-5 ) ! 


A 


he3 
(13-n) 


ov.2  5 
3h  c^ 


A 


(19-n) 


1 


(n-5)! 


3h  c 


A 


(19-n) 


n  =  6,7,8,9,10,11,12 


n  =  13, 14, 15, 16, 17,  II 


V  _     ^ 

2  0    AB 


K 


21    A(A-B) 


K 


22    B(B-A) 


K 


23 


he3 
A  BH 


Kpn  =  hc' 


A(2B  -  10A)(B-A)  +  (15A2-6BA+  B2)  (B  -3A) 


A4(B-A)7 


41 


Koc:  =  he 

25 


15A2  -  6AB  +  B2 


K 


K 


he 


26 


27 


A3(B  -  A)6 
B  -  5A 


he 


A  (B  -  A) 

1 
A(B  -  A)1* 


In  coefficients  Kpj.  to  K„7,  if  A  is  substituted  for  B  and 
B  is  substituted  for  A,  then: 


K   =  K   „ 
n    n-4 


n  =  28,  29,  30,  31 


ho  '   ™ 


c_ 

AB 


K 


K 


K 


K 


K 


K 


41 


43 


44 


45 


46 


47 


D  c  +  jDA1c 

4      p  ? 
2D   +  2AJD 

he3 

"aV 

he3 


32D 


7 


F1F2  +  F^  +  j(P3F2  -  FXF4) 


he 


128D6 
he3 


2     2 
P5Fg  +  P  Fg  +  j(P7Fg  -  F  Pg)- 


2     2 
P6  +  F8 


32D- 


hc 


F9   J'F10 


F 


11 


96D' 


~~ 2     2 
LA1  +  D 


where 


A   «  £ 

Al    2 
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P.    =    26AXD2    -    8A3   +    3WqA. 


F„    = 


i|A1D3    -    4A3D 


F-   =   -    26A2D   +    10D3   +    6cj2D 


aJ   -   6A2D2   +    D 


F   =  32A2  -  32D2  -  12u>2 


F6  =  A3  -  3AlD2 


F„  =  56A1D 


Fg  =  3A2D  -  D3 


Fg  =  -  4A2D  +  3D(A2  -  D2) 


F1Q  =  6AXD2  +  2A3  -  2AXD2 


p  p      p     p  P 

Fll  =  4AiD   +  (Ai  "  D  > 


The  coefficients  K^p,  K^,  K^,  K,-n  and  K,_,  are  respectively 
the  complex  conjugates  of  K^,  s  Ky,  Kiic:»  Kii£  an<3  K^.  „ . 


M]_  =  Re(K4l) 


H 


2  =  Im(Kijl) 


M 


n+1  ~  Re(K^3+1-2n) 


n  =  2,    4,  6, 


43 


M 


=  Im(K I, «  .,  \ 
n       42+^n) 


n  =  4,  6,  8,  10 


2cR 


M 


M 


M 


M 


1 


11    R2  +  R2 
-2cR„ 


12     2     2 
^    R^  +  R^ 

l6E1Rij 


13  uCR^  +  R^) 

i6e1r 

14  w(R^  +  R^) 


N„  = 


coo 


i   R5R6 


N„  = 


coo 
R7R8 


N. 


N,,  = 


"E1(R5R6R18  "  R22R23} 

3R5R6 
"E1^R11R5R6  ~  R9R12) 


N^  =  - 


E1R9 

44 


N^  = 


6R^r5 
5  6 


Nr,   = 


No  = 


N„  = 


_"E1(R7R8R28  " 


R32R33) 


47 
3R7R8 

•E1(R13R7R8  "  R10R14) 

R7R8 

_E1R10 


R7R8 


44 


N 


10   =    6rJh* 


where 


E.    -   -  "I  hc3to 


FL    =    -2oj(A   +   B) 


R0    =    -2(AB    -    oo    ) 


4  4  2    2 

R_    =    R„    +    R_    —    6R„R-| 


R4    =    4R2R^    -    4R^R 


2  2 

Rc    =    A      +    oo 
5 


R6    =   B    -   A 


2  2 

R„    =    B      +    co 


Rg   =    A   -    B 


Rg    =    3A2    -   AB   +    2oo2 


Rio  =  3B?  -  AB  +  2o°2 


Rll    =    B    ~    6A 


R12    ■    9A2    -    4AB   +    5oo2 


R        =    A   -    6B 

Rli4    ■    9B2    -    4AB    +    5oo2 


45 


Rnc:  =  54oo2  +  6B2 
15 

Rl6  =  84A3  -  42BA2 


R.„   =   Rn  CA   -   6B00' 
17  15 


Rl8   "   R16   +   R17 


Rin    =    -    21    A11    +    14    BA3 


R20   =    -   A2(27w2    +    3B2) 

?  4  2    2 

R21    ■    6BAoj      -    lOoo      +    B    a) 


R 


22 


R19    +    R20    +    R21 


R 


23 


12A2    -   6AB   +   6o)2 


R 


25 


54co2    +   6A2 


R 


R 


26 


27 


84B3   -    42AB2 

RocB   -   6Aco2 
25 


R, 


R26    +    R2? 


R 


29 


-      21B11    +    14AB3 


R 


30 


-    (27w2    +    3A2)B2 


R 


31 


6AB032    -    lOoo      +   A2oo2 


R32   -   R29   +   R30   +   R31 
R~_   =    12B2    -   6AB   +   6oo2 


46 


In  coefficients  fi,.,  M  -,  j,  3  N_>  to  N,  n ,  if  -E,  is  substituted 
for  E,  and  3w  is  substituted  for  go,  the  coefficients  M,,, 
M,g3  N, .,  to  N-.  o  can  be  obtained. 


M 


N 


2Aco 


17 


R 


34 


18 


(A2  -  u)2) 


R 


34 


M 


M 


19 


20 


~2E1R38 

2    2 
R38+R37 

—  2E-.  R-oy 

~~ 2      2 
R38  +  R37 


M 


2Aojc 


19 


r*2  4.   2s2 

(A   +  oj  ) 


N 


toe 


20 


2     2 
A^  +  u> 


N21  -^ 


1 

21 


r(A2  +  o)2)Rij0  +  14AR39" 


772  7  278" 

(A   +  oj  ) 


N 


E1R39 


22   3(A2  +  u2)7 


M 


EIR4l 
23  '   (A2  +  co2)6 


N 


N 


E1R42 


24    ^a2     2^5 
6  (A   +  03  ) 


25 


EX(A2  -  Aw2) 
6(A2  +  co2)4 


N 


26 


1 
120 


3A2  -w2 


L 


(A2  +  u)2)3 


47 


AE. 


N 


27         36~0 


(  a2  2x2 

(A      +    co    ) 


N 


28        5040(A2    +   co2 


where 


El   = 


3  ..3 


IT 


he   co 


R 


R 


R 


R 


34 


35 


36 


37 


4A2co2   +    (A2   -   co2)2 


a11    X       4  £A2    2 

A      +    co      -    6A    co 


4Aco3    -    4A3co 


-    2a)(R25    -   R26) 


R. 


R 


R 


R 


39 


10 


41 


-    4coRoc-R0^ 
35    36 

6  ?   4  4   p  fi 

21AU    -   105co   A      +    63co   A      -    3co 


=    -    126A5    +    420co2A3    -    126coV 


S  2    3  4 

3AD    -    10co   A"3    +    3o)   A 


4  2    2  4 

R42    =    5A      -    lOco    A      +    co 


In    coefficients   Min,    M„n,    N        to   NoQa    if   -E,is    substituted 

ly5   20  J   21      2°        1 

for  E1  and  3w  is  substituted  for  co ,  the  coefficients  IVL,  , 
M^pj  N?q  to  N-/-  can  be  obtained. 
2cL. 


M 


13 


23 


2      2 
L13  +  L14 


48 


-2cL 


M 


H 


14 


24  2 

T     +  T 

L13    14 

16E1L18 

25  /T2   .  T2  v 
co(L1g  +  L17) 


M 


26 


16'E -!»--„ 


W(L?8  +  L?7} 


2cojL 


N 


15 


37 


2      2 
L15  +  L16 


N 


N 


N 


38 


39 


40 


N 


41 


-2ctoL,  r 
16 

2      2 

L15  +  L16 

4L6lL63 

L62 

2L32 

L31 

L20  r_  ,., 

'22 


L4L21  +  3L2L3} 


N 


N 


42 


43 


L19L1 


3(L2  +  l|) 
4L61L64 


'62 


2L 


N 


N 


44 


45 


33 


L 


31 
L 


20 


L 


22 


(L3L21  +  3L2L4) 


M 


46 


'19L2 


3(1^  +  h\) 


49 


where 


A   =  2. 
Al    2 


E1  =  -  ^  hcJa) 


2     2     2 

Ll  =  Al  "  D     W 


L2  =  2DA1 


2     2 
L3  "  Ll  "  L2 


L4  "  2LXL2 


5    13    2  4 


L6  =  LiL4+  L2L3 


L7  "  L1L5  "  L2L6 


L8  "  L1L6  +  L2L5 


2     2 
L9  -  4   -   D 


L10  =  A1L9  "  DL2 


Lll  "  A1L2  +  DL9 


2     2 
L12  =  Al  +  D 


L13  =  ~  ^A1(J° 


Ll4  =  ~  2(-L12  -  °°  ^ 


L, c  =  2DL0 


50 


Llg    =    -2DL 


1 


L 


17 


4     4      2   2 
Ll4  +  L13  "  6Ll4L13 


18 


4LmL13  -  4Li4Ll3 


'19 


L 


16D 
L19 


20     D 


21 


3L„  -  L12  +  2oo' 


'22 


2     2 
L4  +  L^ 


'23 


64Df 


'24 


9L9  -  4L12  +  5w 


L 


25 


-  5A1L1  +  7DL2 


'26 


-  7DLX  +  5AXL2 


'27 


L21L24  "  27L2 


'28 


3L2L24  +  9L2L21 


'29 


L23(2DL26  -  L27) 


'30 


"  L23(2DL25  +  L28} 


'31 


2     2 
L5  +  L6 


'32  ~  L5L29  +  L6L30 


L33  ~  L5L30  "  L6L29 
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'34 


84L1q  -  42(A1L   +  DL2) 


'35 


84L11  -  42(AXL2  -  DL9) 


'36 


6L  +  5W 


'37  =  A1L36  +  6DL2 


L38  =  DL36  -  6AlL2 


L39  -  -  6u  A1 


'40 


6co2D 


'4l 


■^34  +  L37  +  L39 


'42 


L35  +  L38  +  L40 


'43 


A1L10  ~  DL11 


'44 


A1L11  +  DL10 


'45 


AiLio  +  DLn 


L 


46 


A1Ln  -  DL1Q 


L47  =  l4L45  "  21L43 


L 


48 


14L46  -  21L44 


L49  =  27oo   +  3L( 


L 


50 


•*  L9L49  ""  3L2 


L 


51 


"  L2L49  +  3L2L9 


52 


r     2  k  2 

Lrn    =    600   L,  n    -    IOoj      +    (o  Ln 
52  12  9 


Lr-«     =     -     0)    Ln 

53  2 


L5lj   "   Li|7   +   L50   +   L   52 


L55   ~   L^8   +   L51   +   L53 


L56  =  6u>  -  6L12  +  12L9 


L57  =  2D(L1Lijl  -  L2L42) 


L58  ~  2D{'L^L^   +   L-^L^) 


L59  "  L5i|L56  "  12L2L55 


L60  "  L55L56  +  12L2L54 


-E. 


L 


61 


768D7 


2     2 
L62  =  Lg  +  L? 


L63  ~  "  L8L58  +  L8L59  "  L7L57  "  L7L60 


L64  ~  L8L57  +  L8L60  ~  L7L58  +  L7L59 


In  coefficients  Mol_,  Moir,  N_„  to  N,,r,  if  -E,  is  subsitituted 

25    26 '   39     ^6 '      1 

for  E-,  and  3^  is  substituted  for  oo ,  the  coefficients   M27> 
Mpo>  Nn„  to  Nf-i,  can  be  obtained. 
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